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Energy and momentum relaxation of heavy fermion in dense and warm plasma 

Sreemoyee Sarkacl and Abhee K. Dutt-MazumdeitU 
High Energy Physics Division, Saha Institute of Nuclear Physics, 1/AF Bidhannagar, Kolkata-700 064, INDIA 

We determine the drag and the momentum diffusion coefficients of heavy fermion in dense plasma. 
It is seen that in degenerate matter drag coefficient at the leading order mediated by transverse 
photon is proportional to (E — fi) while for the longitudinal exchange this goes as {E — n) . We 
also calculate the longitudinal diffusion coefficient to obtain the Einstein relation in a relativistic 
degenerate plasma. Finally, finite temperature corrections are included both for the drag and the 
diffusion coefficients. 



I. INTRODUCTION 



Recent years have witnessed significant progress in understanding the properties of hot and/or dense relativistic 
plasma [l|, |2|. Such studies draw their motivations both from the theory and the experiments. In particular, the 
possibility of creating high temperature quark gluon plasma (QGP) by colliding heavy ions in the laboratory mimicking 
the conditions of microsecond old universe has been a matter of intense research activities in the past decades. Further 

i-G ' impetus to these studies comes from astrophysics where it is important to know the properties of such plasma at high 
density, which, for example, might exist in the core of neutron stars or in white dwarfs. 
fS i' One of the interesting quantities which has assumed special interest recently is the study of partonic energy loss 

D I in relativistic plasma. Several calculations [j-Q have been performed over the last decades to estimate such energy 

tG ^ • loss in a plasma. Similarly, there exists several studies in which momentum diffusion coefficient of heavy fermion has 
been estimated |8l4ll| . These two quantities are of utmost importance to understand the equilibration of fermions 

(^ . in a plasma. So far, these calculations were largely confined to the case of hot plasma with zero chemical potential 
J> ' due to their relevance to the experiments performed at the Relativistic Heavy Ion Collider (RHIC) or the ones to be 
performed at the Large Hadron Collider (LHC). 

'^ ■ There still exists another domain of Quantum Chromodynamic (QCD) phase diagram where the chemical potential 

(/i) might be higher compared to the temperature (T) . This is the re gion of interest of the upcoming experiments on 
compressed baryonic matter (CBM) to be performed at FAIR/GSI |12h14| . Partially motivated by these proposed 

'^ I experiments and partly by another theoretical work on fermion damping rate [l5| , we calculate here the drag (77) and 
the longitudinal momentum diffusion coefficient (S) of a heavy fermion in Quantum Electrodynamic (QED) plasma. 
It is known that the former and the latter are related to the energy loss and the momentum relaxation of the fermion 
in a plasma. Moreover, in equilibrating plasma, these two quantities viz. rj and B are related to each other via Einstein 
relation (ER) which at finite temperature reads as B — 2ETrj. As indicated above, such calculations, for dense (T = 0) 
and/or warm (T <C ^j) plasma are rather limited. In fact, we are aware of only one calculation of energy loss where 
the effects of finite chemical potential has been considered, although the temperature considered there is still high 
C^ I |16| . We on the contrary first consider the extreme case of zero temperature and then incorporate finite temperature 
corrections to our result both for the drag (energy loss) and the diffusion coefficient in the limit fi :$> T. We also 
determine the relationship between 77 and B i.e. ER at zero temperature, which shows some interesting behavior due 
to finite density plasma effect. 

Before we proceed further, it would be worthwhile to draw our attention to [l5|. This is an interesting work in 
many ways. First, it is known that the fermion darnping rate (7) in hot plasma is plagued with divergences which 
cannot be removed by the ordinary screening effect [l7| . This is because, the magnetic interaction is screened only 
dynamically [13 and the problem remains for the static photons (or gluons in QCD) . Therefore, to obtain finite result, 
a suitable resummation has to be performed. This was first done in [l^ [201 ■ -R-^^- UM shows that at zero temperature 
due to Pauli Blocking, finite result can be obtained without performing further resummation. This is consistent with 
the conclusion drawn in |2l| . Secondly, in the relativistic plasma 7 is dominated by the magnetic exchange and is 
proportional to {E — /x), while the electric photon exchange gives a contribution proportional to [E — /x)^. Here, it is 
important to note that the dynamical screening in the transverse sector enhances the damping rate compared to its 
longitudinal counterpart. It might be recalled also that for non-rclativistic Coulomb plasma the damping rate goes 
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as {E — ^Y [22[- Thus, it would be interesting to see how do the drag and the diffusion coefRcient depend on [E — /i) 
in degenerate plasma. 

It is known that at finite ternperature the calculation for the energy loss and diffusion coefficients are plagued 
with infrared (IR) divergences [l4|. To deal with this problem, in hot plasma one separates the integration into 
two domains: one involving the exchange of hard photons (or gluons) i.e. the momentum transfer [q) ~ T and the 
other involving soft photons (or gluons) when q '^ eT (e <^ 1). In case of the former, one uses bare propagator 
and introduces an arbitrary cut off (q*) [23| parameter to regularize the integration. For the latter, on the other 
hand one uses the Hard Thermal Loop (HTL) corrected propagator. These two parts, upon addition, yield results 
independent of this intermediate scale. In case of degenerate plasma, one also encounters similar infrared divergences 
and following ref. [l5| one can proceed along the same way as finite temperature (using Hard Dense Loop (HDL) 
corrected propagator) and show that both for the drag and diffusion coefficient the final result becomes independent 
of the intermediate cut-off parameter. This however, as we shall see, is not required in case of dense plasma. Here, 
the dominant or the leading order contribution comes entirely from the soft sector and the hard photon exchange 
contributes only to the higher order. It might be mentioned here that although we calculate these quantities for QED, 
with appropriate color factors the results can easily be extended to the case of QCD with the addition of one more 
diagram involving triple gluon vertex [J]. Furthermore, it might be noted that the quark energy loss calculations in 
general should also include Bremsstrahlung radiation of the gluons. However, in the present context we are concerned 
with only the two body scatterings and therefore restrict ourselves to the coUisional energy loss alone. 

Furthermore, expressions derived for the degenerate plasma, wherever possible, have been directly compared with 
their finite temperature counterparts (with zero chemical potential). This brings the similarities and the differences 
of these two extreme scenarios into clearer relief. 

The plan of the paper is as follows. First in section II we calculate drag and diffusion coefficients in degenerate 
plasma and discuss about Einstein relation. In section III, the finite temperature corrections have been incorporated 
both for 77 and B. The results are then summarized in section IV. An appendix has also been added to understand 
the origin of difference in ER in cold medium than from that of hot plasma. 

II. HEAVY FERMION AT ZERO TEMPERATURE 

A. Drag coefHcient 




FIG. 1: Fermion self-energy with resummed photon propagator. 

In this section we first calculate the drag coefficient of a heavy fermion in a degenerate QED plasma. For this we 
consider scattering of a heavy fermion having energy (E) (which we assume to be hard), with the constituents of the 
plasma viz. the electrons. Incidentally, this drag co-efficient (77) is related to the to the energy loss by the following 
equation: 

Evi V dx 

where, v^ = -^ is the velocity of the incident fermion, {dE /dx) is the energy loss and p is the three momentum of 
the incident fermion. Thus, the calculation of the drag coefficient boils down to the calculation of collisional energy 
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loss in a plasma 0-13 • Now the energy loss (dE/dx) can be obtained by averaging over the interaction rate times the 
energy transfer per scattering uj and dividing by the velocity of the incoming particle [3| , 



da; Vj 



(2) 



This expression is quite general and valid for both the finite temperature and/or density where only the phase space 
will be different due to the modifications of the distribution functions depending upon the values of ^ and T. 

The scattering rate, which is essential for the calculation of r] as evident from Eqs.(l2]) and ^ is related to the 
imaginary part of the fermion self energy (E) by the following equation |24| : 



r(S) = -^Tr [ImI](po + "/ ,p)(f + m)] 



Po=E 



(3) 



In the last equation, m is the mass of the incoming heavy fermion. The full fermion self-energy represented in Fig.(IT]) 
can be written explicitly as: 



E(P) = e^T^ y 7^'^'^ Sfiiiun - ujs),p - q)7. A^,(^L^„ q) 



(4) 



where, po = iuj„ -\- n, q^ — itOs- a;„ = 7r(2n -|- 1)T and Wg = 2'ksT are the Matsubara frequencies for fermion and 
boson respectively with integers n and s. After performing the sum over Matsubara frequency in Eq.Q, iujn + /i is 
analytically continued to the Minkowski space icj„ + fi ^ po + irj, with 77 — >■ 0. The blob in Fig.([T]) here represent 
HTL/HDL corrected photon propagator which is in the Coulomb gauge is given by [l| , 



with, Pl^ = {5,j - m,), <r = qV|q|, Ph = P^ = ^oo = O and A,, At are given by Q, 

-1 



At(<?o,'?) = 



cf + Ui 



-1 



qt> -q^-^t 



For subsequent calculations it is convenient here to introduce the spectral functions pi^t \^'- 

Pi.t{qo,(l) = 2ImAz,t(go-t-i?/,q)- 



(5) 

(6) 
(7) 

(8) 



At the leading order these are derived from the one-loop photon self-energy where the loop momenta are assumed to 
be hard in comparison to the photon momentum [25l . l2a | . In the literature the formalism is known as the HTL/HDL 
approximation as discussed in |l|. 



pi{qo,q) 



Pt{qa,q) = 



27rTO|,a;e(l-a;2) 



92 + m]j ( 1 - f In 



x+l 
x-1 



27rm|, v'jx{l- x'^)Q{l - x^) 



2q'^{x'^v'l — 1) — m?jjx'^v'l I 1 



(i-^-^) 



In 



x+l 
x-1 



rn'^v^ 7r^ x'^ (1 — x-^ )^ 
4 



(9) 



2x 

where, Vf is the Fermi velocity and x ~ -22-. For a ultrarelativistic plasma {vf — > 1) the Debye mass is mjj 



1^\ 



In Eq. ^ , fermion propagator has the following spectral representation with the notation k = (p — q) [l| , 



S'/(iw„,k) 



dfco i^ + m)pfiK) 



27r kc) — 



- lOJn - M 



(10) 



Hence, for Pf{K) we use the free spectral density given by, 

Pf{K) = -^m^ - Ek) - S{ko + E,)] . (11) 

One can take the imaginary part of the Eq.Q to calculate the scattering rate with the help of Eq.Q. For the 
calculation of the drag coefficient, one then inserts the energy exchange to in the expression of F and calculate dE/dx 
from Eq.(l2]) to obtain, 

dE ne^ f d'q f°° dfco ., ^ p dqo 

X (l + n(go)-n(fco))(5(£'-fco -go) 

X [poA;o + P •k + m^]p(((7o,'?) 

+ 2[pofco-(p-q)(k-q)-™^](0t(go,g)- (12) 

It is to be mentioned here that the scattering process involves space like photons. Hence, here only the cut of the 
spectral function contributes. In the above equation n and n are the Bose-Einstein and the Fermi-Dirac distribution 
functions: 

From now onwards in this section, we exclusively focus on the ultradegcnerate plasma. The finite temperature 
corrections which might be important for dense and warm plasma will be incorporated in the next section. For T = 0, 
/i 7^ limit, (1 + n{qa)) = Q{qo) and n(ko) — Q{fi ^ E + qo), where O represent the step function. These functions, as 
we shall see, restrict the phase space of the qo integration severely. The zero temperature spectral functions pu now 
involve Debye mass mjj = e'^VffP/n'^. 

Note that, the delta function in Ea.(|12p sets qo = qvcosO and the theta functions impose further restrictions on go- 
We consider quasiparticles, with velocity close to Fermi velocity, which undergo collisions with the particles near the 
Fermi surface. Hence, we can make an approximation here as u sa wy. 

Now, consider the case of hard photon exchange where the medium effects on the photon propagator can be ignored. 
In this case using the bare propagator we get, 

/ dE\ e^mjy 




e'^iE - pfmjy 1 

2ATTVf Vf 



''i- + ^)l^- (14) 



This actually is the leading hard contribution that comes from the diagram, when, the blob of Fig.([T|) is replaced 
with one fermion loop. Evidently, the above integral is infrared divergent and unlike the finite temperature here 
higher powers of q appear in the denominator. We shall remark on this later once we have expressions both for 77 and 
B. 

To deal with this infrared divergence in the soft domain, one uses HDL corrected photon propagator [25l l26j given 
by Eq.Q, with Debye mass mjj = e^Vjjj? /'k'^ as mentioned earlier, to obtain, 

/ dE\ e? /■ d^a 

[-^)ijE) . ^y (^,./cos.(e(,o)-e(,-i. + ,o)) 

X Q{q* - q){pi{qQ, q) + v){l ~ cos^ e)pt{qo, q)} 

dqq I dqaqa{pi{qa,q) + {vj - -^)pt{qo,q)} ■ (15) 



Stt^vj J jj J q 

The integration domain(£') above is limited by the Q functions. 



D : <qo< E - p; 

qo <q<q*- (16) 



With these we get, 



dx ) 



soft 



{E) 




(17) 



where, 



n I \ 1 X l + x 



Qt{x) ^-Qi{x) + 



l-x^ 



(18) 



We are mainly interested in the energy loss of a quasiparticlc which is present close to the Fermi surface, hence, 
(£^ — ^) <C rriD is the physically interesting region where the quasiparticlc concept is meaningful. The denominator 
of the Ea. (fT7|) can now be expanded in powers of q^. We replace s* = {q* /mo)'^ and compute separately the electric 
(Z) and the magnetic part (<), 



dE\ 



V dx 

/ dE 

V dx 



e^E-fif 



ds 



oft iSmriDV^j: Jq y/s{s + 1)^ 

dqodq- 



t 

soft 



2 ? ? 
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(19) 

(20) 



After explicit calculation, the electric and magnetic contributions to the expression of energy loss take the following 
form. 



dE\ 
dx / 



e^iE-^if e'^{E- fj.)^mj 



soft 



96w?7Tl£) 



12'Kv].q*^ 



(21) 



dE\ 
dx ' 



e^{E-iif e^{E-nfvf7 



'n 



soft 



487rw/ 



lAAirq 



*3 



(22) 



It is worthwhile to note here that the leading order terms in the last two equations arc finite and independent of 
the cut-off parameter. Here the q* dependent term appear only at O(e^). Therefore, we write. 



dE\ I 
dx ) 



soft 



e^jE^fi)- 
96u?m£) 



+ Oie% 



(23) 



dE\ 
dx I 



soft 



A&TTVf 



+ 0{e'). 



(24) 



So far we have discussed about the soft part and have seen that in the limit q* — >■ oo the cutoff parameter dependent 
term trivially vanishes. Similarly, if we recall the expression for the hard part i.e. Ea. (jl4p . after performing the 
integration in the limit [g*, ^\ we get. 



dE\ 
dx J 



e'iE-nYmjj 1 vj 



hard 



727TVf 



(7^ + ^) 



^*3 



0{e\ 



(25) 



Clearly it fails to contribute at the leading order where the entire contribution comes from the soft sector. This 
is a distinctive feature of degenerate plasma not encountered at finite temperature (/x = 0). There both the hard 



and the soft part contribute to the leading order in e^ and the divergence is only logarithmic. To deal with such 
divergences in hot plasma one invokes Braaten and Yuan's prescription ^23j| where an intermediate cutoff is introduced 
to separate the hard and the soft domains. It is seen that such an intermediate cutoff parameter disappears from the 
final expressions when both the contributions are added. At zero temperature, a similar approach was adopted for 
the calculation of fcrmion damping rate [la] where it was shown that such cancellation takes place also in degenerate 
plasma. It is obvious from Eqs. ((2T1 [22l [25]) that same thing happens for 77 also. 

From Eq. (f25|) it is clear that the result obtained from the hard region is suppressed with respect to the soft one 
fEg. ipT]) and ([H])). Hence, the whole contribution to leading order comes from the soft sector alone. The final 
expression for drag coefficient at zero temperature becomes. 

The first term above corresponds to the electric photon and the latter to the magnetic one i. e. I or t mode behaves 
differently. The dominant contribution to rj comes from the magnetic sector in the ultrarelativistic case v/ — >■ 1 and 
the electric sector when Vf -^ I. Results for the light fermion can be obtained from Ea.p6| with the substitution of 
Vf -^ 1. 

B. Diffusion coefficient 

Apart from rj, the quantity which could be of imp ortance in the study of heavy fermion propagating in the plasma 
is the momentum diffusion coefficient (i?y) |8l4llj. In fact, we know for Coulomb plasma rj and the longitudinal 
momentum diffusion coefficient (B) are related via ER. Momentum diffusion coefficient Bij can be defined as follows 

B,,^fdrq,q,. (27) 

Decomposing Bij into longitudinal {Bi) and transverse components (Bt) we get the following expression, 

B.,,^BS^J-^)+Bl^. (28) 

These coefficients Bit are the longitudinal, transverse squared momentum acquired by the particle through collision 
with the plasma. Using the above definition, like the energy loss (Eq. (jl2p ). longitudinal momentum diffusion coefficient 
{Bi = B, suppressing the index I) can be written as follows, 

X {I + n{qn) - n(ko))5{E - kn - qn) 

X [pofco + P ■k + m^]pi(qo,'7) 

+ 2[pofco-(p-q)(k-q)-m2]pt(go,g). (29) 

Here, qu = qcos9. For the exchange of hard photons using the bare propagator we obtain: 



s.^^±4^i.,r.4 



^^•"^<^-<(^4)/^^ (30) 



32'Kvj vj 2' J q^' 

Comparing Eqs. (jl4p with pop it is seen that like ry, B is also infrared divergent involving fourth power of q in the 
denominator. At finite T case, both the quantities are proportional to dq/q at the leading order [1^. To understand 
the origin of this difference, we focus on the go integration. It is shown in the appendix that in medium at finite T, 
there involves quadratic power of go in both cases with the limits —vq to +vq giving rise to a term proportional to q^ 
in the numerator. This cancels with some of the powers of q coming from the propagator. Whereas in cold matter, 
from Eqs. (|14p . pop we find that the same integrations appear with q^ and ^q in the numerator while the limits are 
independent of q forbidding the cancellation with g's coming from the propagator as before. We note here that the 



drag and diffusion coefScicnts are related through B = — 4^ fj when we deal with the bare propagator. We shall 
see in the next paragraph that same powers of go appear in the numerator when one takes the plasma effects into 
account but such common scaling behavior is lost. 

The infrared divergence of Eq. (pO| can be removed by using the dressed photon propagator [2^, [2g and providing 
the upper cut off as in the case of 77. With the HDL corrected propagator one gets, 



B 



soft 



(E) 



2v} 



(2^)3 



g2(e(go) ^ei^i-E + qo)) 



X e{q* - q){pi{qo, q) + vj{l - cos^ 0)pt{qo, q)} 
dqq ' 'J"-"2r„./„_ „^ , /„.2 



87r2f2 j^ 



ll 



dqoqo{pi{qo, q) + {vf j)pt('7o, g)} 



(31) 



here the integration domain (D) is same as before. So, 



B\sott{E) 




2q^+mlv}Qt{f) 



1 2 



2q2 



vWq^ /■ 



4g2 



(32) 



Since, we know from energy loss that dominant contribution to the expression comes from the soft region alone we 
write the expression for B as follows. 



B 



e^{E~ij)^ ^ e^jE-pf 
128TODWf 



727ri;2. 



+ 0{e^), 



(33) 



which is finite. Now from Eas. (j26|) and (|33p it can be seen that there is no common scaling factor between rj and B. 
But as ER is formulated in the region where u^: ^ 1, in this nonrelativistic region exchange of the magnetic photons 
are suppressed in comparison with the electric one. Hence, considering only the electric part we get the same ER, 
B = — \^ 11 as in the case of bare perturbation theory. 



III. FINITE TEMPERATURE CORRECTION 



The results of the previous section can easily be extended to the case of a hot and dense (T <C /x) plasma. This could 
be relevant for heavy ion collision to be performed at GSI where the chemical potential is expected to be much higher 
than the temperature. Now, while calculating the soft part we replace the zero temperature distribution functions 
with the finite temperature one in Eq. (|15p and write. 



dx J 



e^ [ d^q 

A^) - t;— 77rT3<l^f(^os9{l + n(qo)-n{E-qo)) 

soft ZVf J [ZTT) 

X e(<7* - q){pi{qa, q) + v}{l - cos^ e)pt{qa, q)} 



With small T and large p, the above equation can be calculated according to Ref.|2'/ 
write any function g(e) along with the fermion distribution function as follows. 



(34) 
In this approach we can 



9ie) 



e/3(^-A') + 1 io 6 



-de 



9\P^)- 



The contributions connng from soft {I and t) using Eas. (j341 135[ 1^ are found to be given by: 



(35) 



dE\ 
dx ) 



e^E-pf e^{E^pfml e^E ~ p)T^Tr^ e^{E - p)T^ 7x1^,11 



soft 



96w?7Tl£| 



72Trq*3vj 



96z;?TOi3 



72q*Hj 



(36) 



dE\t e^iE-nY e\E- ^yvfmjj e^nT^ e^{E ~ fi)mjyVfT^7r 



dx I 



soft A8nvf 1447rq*3 72vf 144q 



*3 



(37) 



We see from the above two equations that the terms containing separation scale are subleading in comparison with 
the others, the same behavior obtained in the zero temperature case also. The term with the bare propagator comes 
as, 



dE\ 
dx J 



hard 12nVf u? 2' 



1 1 



/ 
^ (^ + -)—--• (38) 



1 1 

q*3 ^3 



The above expression is also suppressed in contrast to the soft one. Hence, with finite temperature correction, drag 
and diffusion coefficients become, 

^ e'{E-^,r . e'jE-^,)-' e^Tr^T^jE - ^,) e'nT' , 

^~ QGniDvjE '^ ASnvjE 96mDvjE 72v)E^ ^ '' ^ ' 

128mDvj 72TTVJ 6AmDvj A8v} + ^^^ ^ ^^"^ 

One notes here with the thermal correction the ER cannot be established even for the electric sector alone. 

IV. SUMMARY 

In this work we calculate the energy loss and momentum diffusion of heavy fcrmion in dense and warm QED matter 
and highlight some of the differences that exist between the hot {jj. = 0) and the cold {T = 0) plasma. Unlike finite 
temperature, where one encounters logarithmic divergences in calculating t] or B, here we come across non-logarithmic 
divergences. Furthermore, we see that at the leading order in coupling, the entire contribution comes from the soft 
sector and this is finite i.e. the physics here is dominated by the excitations near the Fermi surface. The exchange of 
hard photons on the other hand contribute only at O(e^). It is to be noted that in a thermal medium with vanishing 
chemical potential both the soft and hard photons or gluons for QCD matter contribute at the same order. Moreover, 
for ultrarelativistic particles both ry and B receive dominant contributions from the magnetic sector while the electric 
parts are found to be subleading in {E — fi). This is consistent with fermion damping rate calculation [1^ [2l|, ^^ 
in degenerate plasma. Quantitatively, we find that for the transverse or magnetic interaction rj is proportional to 
{E — jj,)^ while for the electric interaction, it goes as {E — n)^. Similar differences for B is also seen where one more 
extra power of {E — fj.) involved in each case. The other important finding of the present investigation is the ER for 
the drag and diffusion coefficient. In hot plasma, it is known that B ~ 2TE'q [9l4ll|. At zero temperature, we find, 

B = — A~'^ 11 by considering only the bare propagator i.e. when we do not take the plasma effects into account. 
However, we see that this common scale behavior is lost for soft photon exchange where the plasma effects are included 
and both the magnetic and electric contributions are retained. However, by retaining only the electric contribution 
for the cold plasma, one arrives at the same relations as obtained by using the bare propagator. For T ^ fi again we 
see that, with plasma effects incorporated, rj and B fail to show such common scale behavior even when the magnetic 
interaction is ignored. 

As a last remark, we note that here the entire calculation has been done for QED plasma. This can easily be 
extended to QCD matter with appropriate modifications like inclusion of diagrams involving three gluon interaction 
and proper vertex factors coming from the QCD color algebra. Such studies are in progress and shall be reported in 
future. 
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VI. APPENDIX 

To understand the difference of tlie results between tlie cold and hot plasma, we first recall the expression for the 
drag coefficient (77): 

Evi \ ax J 
The above relation with Egs. pU]) . ([TT|) and p^ can be further simplified to yield, 

'^ - Q 2 217 / '^gg / '^<ioqo{l + n{qo)-n{ko))e{q* -q){pi{qo,q) + {v'^ - -^)pt{qo,q)}- (41) 

OTT V -tj J J-vq Q 

The corresponding expression for the diffusion coefficient from Ea. (P5)) is, 

^ - ^^ f'^11 I dgogo'(l + «(9o) - n(fco))e(g* - q){pi{qo.q) + («' - %)pt(go,g)}- (42) 

^'' ^ J J — vq q 

In the high temperature limit, (1 + n(go) — n{E — go)) — ~ + 5- It is to be noted that the above integration 
limits are symmetric in go- Hence, for the drag, the factor of i and for the diffusion — contribute. Inserting these in 
Eqs.dm and ^ we get, 

g2 r rvg „2 

^ - 1^ 2 277 / dgg / dqoqoOiq* - q){pi{qo,q) + {v'^ - -^)ptiqo,q)}, (43) 

iOTT V hi J J -vq q 

and 

g22i /• ^vq 2 

g - o 2 2 / ^^gg / dqo(?o6(g* -g){Pi(go,'?) + (w^ - -|)pt(go,g)}- (44) 

OTT U J J —vq q 

In case of bare interaction, one can show that both ?/ and S is proportional to dq/g [lO[ and even without performing 
the integration B = 2TErj. If we compare the Egs. pS)) and (|33|) with Eqs.([T5|) and ([31]) we find that the q^ integration 
for cold matter is not symmetric, and the limits are independent of q. Here, lies the difference of cold and hot plasma. 
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